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ABSTRACT 


In  this  report  we  describe  optical  disks  that  store  data  holographically  in  three  di¬ 
mensions  using  either  angle  multiplexing  or  wavelength  multiplexing.  Data  is  stored  and 
retrieved  in  parallel  blocks  or  pages,  each  page  consisting  of  approximately  one  million 
bits.  The  storage  capacity  of  such  disks  is  derived  as  a  function  of  disk  thickness,  pixel 
size,  page  size,  and  scanning  parameters.  The  optimum  storage  density  is  approximately 
120  bits/pm2. 
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1.  Introduction 


Holographic  storage  of  data  in  3-D  media  such  as  photorefractive  crystals  can  provide  high  storage 
density  and  fast  parallel  access  to  the  stored  information.  Such  memories  were  investigated  very  extensively 
in  the  early  60’s  [1-5].  Even  though  these  early  efforts  produced  remarkable  results  they  never  found  practical 
application  largely  because  of  material  limitations  such  as  low  sensitivity,  fanning,  and  hologram  decay. 
Interest  in  holographic  3-D  memories  has  been  revitalized  in  recent  years  for  a  variety  of  reasons,  including 
significant  advances  in  recording  materials  (photorefractive  and  photopolymers),  dramatic  improvements  in 
all  optical  devices  (lasers,  spatial  light  modulators,  detectors,  etc.),  and  most  significantly,  the  emergence  of 
applications,  such  as  neural  networks,  machine  vision,  and  databases,  that  can  make  use  of  the  capabilities 
of  holographic  3-D  memories.  The  theoretical  upper  limit  on  the  storage  density  is  V /A3,  where  V  is  the 
volume  of  the  hologram  and  A  is  the  wavelength  of  the  light.  This  limit  is  in  the  order  of  1012  bits  per  cm3, 
however  in  practical  systems  only  109-1010  bits  per  cm3  is  achievable  due  to  the  finite  numerical  aperture  of 
the  optical  system  that  transfers  the  data  into  the  optical  system  and  the  dynamic  range  of  the  crystal.  For 
example,  103  can  be  superimposed  at  the  same  location,  each  hologram  consisting  of  103  x  103  pixels,  giving 
a  total  memory  of  109  bits  per  location.  The  practical  usefulness  of  such  a  memory  must  be  considered 
in  light  of  the  fact  that  electronic  RAM  memory  chips  currently  under  development  have  storage  density 
of  64  million  bits  [6].  Twenty  such  chips  mounted  on  a  single  board  could  match  the  storage  capability  of 
the  volume  hologram.  In  order  to  build  a  mass  storage  medium  that  is  not  threatened  by  semiconductor 
memories,  we  must  construct  holographic  memories  that  have  a  capacity  much  larger  than  10®  bits. 


Figure  1.  The  3-D  holographic  disk  (HD). 


In  this  report,  we  present  a  spatially  multiplexed  3-D  holographic  storage  scheme  that  w»  refer  to  as  the 
3-D  holographic  disk  (HD)  [7],  As  in  all  spatial  multiplexing  schemes,  the  most  crucial  component  of  the 
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system  is  the  scanning  mechanism  that  steers  the  readout  mechanism  to  different  locations  on  the  disk.  In  our 
system,  spatial  multiplexing  is  done  in  a  disk  configuration  with  the  rotation  used  to  access  different  recording 
locations,  as  shown  in  Figure  1.  Two  light  beams  (a  signal  and  a  reference)  interfere  inside  the  photorefractive 
crystal  to  create  a  phase  grating  via  the  photorefractive  effect.  Multiple  holograms  are  recorded  at  the  same 
location  by  changing  the  reference  beam  angle  (angle  multiplexing)  or  by  changing  the  wavelengths  of  the 
reference  and  signal  beams  (wavelength  multiplexing).  Because  of  the  Bragg-matching  requirement  of  volume 
holograms,  individual  holograms  can  be  read  out  by  changing  the  direction  of  the  reference  beam  (for  angle 
multiplexing),  or  the  wavelength  of  the  reference  beam  (for  wavelength  multiplexing).  We  will  assume 
throughout  that  the  image  beam  is  at  normal  incidence  on  the  crystal. 

The  main  result  of  this  report  is  the  derivation  of  the  storage  capacity  of  3-D  HDs  as  limited  by  geomet¬ 
rical  constraints.  We  show  that  a  3-D  disk  of  approximately  1.5  cm  thick  has  storage  density  approximately 
equal  to  120  bits/ /an2.  Thus,  a  3-D  HD  stores  the  equivalent  of  more  than  a  hundred  conventional  2-D  disks 
of  the  same  area. 


2.  Angle  Multiplexed  Holographic  Disk 

In  this  section  we  address  the  following  question:  What  is  the  maximum  number  of  bits,  JV,  that  can 
be  stored  in  a  3-D  HD  of  area  A  using  angle  multiplexing?  We  will  show  that  in  order  to  maximize  N 
we  must  properly  select  the  thickness  of  the  HD  (£),  the  magnification  of  the  optical  system  that  transfers 
the  data  to  the  disk,  and  the  angles  of  incidence  for  the  reference  beam.  In  what  follows  we  derive  these 
optimum  parameters.  The  limits  to  storage  capacity  in  this  report,  are  due  to  geometrical  constraints.  The 
dynamic  range  of  the  recording  material  imposes  a  limit  on  storage  density  independently.  We  will  see  that 
the  capacity  due  to  the  geometric  constraints  is  more  restrictive  than  the  material  limitations  in  the  3-D  HD 
system. 

We  can  express  N  as  follows: 


N  =  N.NeNl.  (1) 

In  the  above  equation  N,  is  the  number  of  separate  locations  on  the  disk  where  holograms  are  superimposed, 
Ng  is  the  number  of  holograms  that  are  angularly  multiplexed  at  the  same  location,  and  jV2  is  the  number 
of  pixels  in  each  stored  hologram.  We  will  derive  an  expression  for  each  of  the  three  quantities  and  then 
maximize  their  product  with  respect  to  the  various  parameters  of  the  system. 


2a.  Maximum  Number  of  Angularly  Multiplexed  Holograms 

We  derive  an  expression  for  the  maximum  number  of  holograms,  Ng ,  that  can  be  angularly  multiplexed 
at  a  single  location.  In  the  following  analyst  data  is  stored  by  recording  either  reflection  or  transmission 
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holograms.  The  reference  beam  is  a  planewave  whose  incident  angle  is  Or.  The  signal  beam  can  be  considered 
as  a  superposition  of  planewaves  that  spans  a  range  of  angles.  To  calculate  Ng  we  must  first  calculate  A Or, 
the  full  width  of  the  angular  selectivity  of  each  hologram,  which  we  take  to  be  the  angular  separation  between 
adjacent  holograms.  An  approximate  expression  for  AOr  is  [8] 


A  Or 


8A  cos  Os 
nwL  |  sin  (Or  +  0S)| ' 


(2) 


where  A  is  the  wavelength,  L  is  the  thickness  of  the  hologram,  n  is  the  index,  and  Os  is  the  incident  angle  of  the 
central  planewave  component  of  the  signal  beam  (see  Figure  2).  For  transmission  holograms  0  <  |0/f|  <  tt/2, 
and  for  reflection  holograms  n/2  <  |$/?|  <  n.  The  signal  beam  is  assumed  to  be  in  the  range  0  <  |  <  v/2. 


Figure  2.  Recording  geometry. 


Eq.  (2)  is  only  an  approximate  estimate  for  the  angular  selectivity  of  the  entire  grating  since  different 
planewave  components  have  different  A  Or.  However,  Eq.  (2)  is  commonly  used  for  setting  the  angular  sep¬ 
aration  between  reference  beam  angles.  The  cross-talk  resulting  when  holograms  are  angularly  multiplexed 
in  this  way  has  been  calculated  recently  [9]. 

To  calculate  the  number  of  holograms  that  can  fit  into  a  range  of  reference  beam  angles  Or  spanning 
from  to  02  (each  hologram  being  separating  from  its  adjacent  holograms  by  a  corresponding  A  Or),  we 
observe  that 


8A 

|sin(0fl  +0s)\A0r  =  ^jcos0s,  (3) 

which  is  valid  for  all  possible  angles  Or.  If  we  add  together  Ng  —  1  such  equations,  one  for  each  value  of  Or, 
and  approximate  the  left  hand  side  of  the  summation  by  an  integral,  we  obtain  the  following  expression: 
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/•«> 

/  |  sin(0fl  +  0s)|d0fl  = 

Jg  i 


8X(Ng  —  1) 


COS  0S- 


Solving  for  the  number  of  reference  angles,  we  get 


Ng  =  1  + 


(nn  L\ 

~&r) 


I  COS(0g  4~  6\  )  —  COs(#5  +  02  )| 
cos  Os 


where  it  is  assumed  that  either  0  <  6$  +  6\  <  6s  +62  <  tr.  or  —n/2  <  6s  +  6\  <  6s  +  62  <  0.  Physically, 
this  means  that  the  reference  beam  is  always  to  one  side  of  the  signal  beam.  The  above  calculations  were 
carried  out  for  angles  inside  the  crystal.  We  can  use  Snell’s  law  to  convert  to  angles  outside  the  crystal. 

In  the  following,  we  will  assume  that  the  image  beam  has  normal  incidence  (6s  =  0).  In  this  case, 
Eq.  (5)  becomes 


Ng  =  !  +  I  cos  *  -  cos02|, 


where  we  have  0  <  6\  <  02  <  Jr/2  for  transmission  holograms.  We  can  increase  Ng  by  a  factor  of  2  by 
recording  a  second  set  of  angularly  multiplexed  holograms  in  the  range  —  6\  to  — 02,  since  the  number  of 
holograms  that  can  be  angularly  multiplexed  in  the  same  range  of  angles  is  equal  to  the  expression  in  Eq.  (6). 
It  is  also  possible  to  simultaneously  angularly  multiplex  reflection  and  transmission  holograms,  as  shown  in 
Figure  3.  Therefore,  the  geometric  limit  on  the  total  number  of  holograms  that  can  be  superimposed  in  the 
same  location  is  4  times  the  expression  in  Eq.  (6)  . 


REFLECTION 

HOLOGRAMS 


TRANSMISSION 

HOLOGRAMS 


SIGNAL  BEAM 


Figure  3.  Angular  multiplexing  by  reflection  and  transmission  holograms  from  both  sides  of  the 
signal  beam. 
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2b.  Spatial  Multiplexing 


The  number  of  non-overlapping  spatial  locations  on  a  disk  with  area  A  is 


N.  =  -=—r  (7) 

a  ww' 

where  a  =  w  x  w'  is  the  area  of  each  location.  To  determine  w  and  w'  we  need  to  take  into  account  the 
fact  that  the  stored  images  can  be  in  exact  focus  at  only  one  plane  in  the  volume  of  the  crystal.  As  the 
thickness  of  the  crystal  increases,  the  area  occupied  by  the  defocused  image  at  the  surface  of  the  hologram 
also  increases.  Moreover,  the  size  of  the  area  that  is  illuminated  by  the  off-axis  reference  beam  increases 
in  one  dimension  as  the  crystal  thickness  and  the  angular  sweep  increase.  We  will  derive  expressions  for  w 
and  w'  with  reference  to  the  geometry  of  Figure  4.  We  assume  that  the  images  to  be  stored  are  at  normal 
incidence  and  are  focused  at  the  middle  of  the  crystal.  We  can  calculate  the  extent  of  the  defocused  image  on 
the  surfaces  by  tracing  the  rays  corresponding  to  the  highest  spatial  frequency  of  the  focused  image.  Let  6  be 
the  resolution  or  pixel  spacing  of  the  focused  image.  Then  the  maximum  spatial  frequency  is  approximately 
1/6,  corresponding  to  a  diffracted  plane  wave  traveling  at  an  angle  9  =  sin-1  (A /n6).  We  use  the  ray  optics 
approximation  to  trace  this  maximum  spatial  frequency  component  and  obtain  the  size  of  the  defocused 
image  at  the  crystal  faces: 


Figure  4.  Angle  multiplexing:  extra  area  taken  up  by  defocusing  and  reference  beam  angle  change. 


w-  N,,6  +  Lt&n9  =  Np6+  .  ==.  (8). 

yj{n6/\y  -  1 

As  shown  in  Figure  4,  in  order  for  the  reference  beam  to  fully  illuminate  the  volume  of  the  crystal  that 
the  signal  beam  occupies,  it  must  illuminate  a  width  larger  than  w  in  the  direction  of  reference  beam  sweep. 
From  the  geometry  of  Figure  4,  this  width  is 
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w'  =  tv  +  L  tan  82  . 


(9). 


The  overall  area  that  must  be  devoted  to  each  recording  location  is  therefore 

a  =  ww'  =  w(w  +  L  tan #2),  (10) 

where  w  is  given  by  Eq.  (8). 

2c.  Optimum  Np>  9 1,  and  6 

We  can  now  write  an  expression  for  N ,  the  total  number  of  bits  stored,  by  using  Eqs.  (7)-(10)  and 
Eq.  (1): 


N  =  AN; 


Nt 


ww 


1  +  ^(cosfl,  —  cos  82) 


=  AN- 


_1A 


NJ i  + 


vWM)5- 


Nvb  +  L  +  L  tan  9-> 

V  v/lnd/AP-l 


(11) 


We  wish  to  maximize  the  above  expression  by  optimally  selecting  Np,  £,  9 1,  8,,  and  6,  which  are  the 
parameters  we  can  control. 

First  of  all,  we  note  that  N  decreases  monotonically  as  9\  increases  (in  our  analysis  0  <  fli  <  x/2), 
therefore  8 1  =  0  is  the  optimum  value.  However,  since  the  angular  selectivity  is  very  poor  around  =  0, 
in  practice  the  minimum  angle  of  the  reference  is  set  at  #i  s;  10°  inside  the  crystal.  Next  we  consider 
the  optimum  number  of  pixels,  Nv.  Taking  the  derivative  of  N  with  respect  to  Np  shows  that  N  is  a 
monotonically  increasing  function  of  Np.  This  result  confirms  our  intuition  since  the  increase  in  the  disk 
area  required  to  store  the  holograms  due  to  defocusing  and  angular  multiplexing  can  be  thought  of  as  an 
“edge”  effect.  The  use  of  larger  images  implies  fewer  recording  locations  on  the  same  disk  area,  and  hence 
fewer  edges.  In  practice  Np  is  limited  by  the  number  of  pixels  of  the  spatial  light  modulator  (SLM)  to 
approximately  Np  =  1,000.  For  the  rest  of  this  section,  we  will  consider  0i  and  Np  as  given  and  fixed. 

The  determination  for  the  three  remaining  variables  (L,  82,  and  6)  is  more  difficult.  We  first  consider 
the  optimum  pixel  size  6.  For  a  given  L,  N  is  maximized  with  respect  to  6  when  w  is  minimized  with  respect 
to  6.  To  find  the  optimum  6  it  is  convenient  to  write  to  as 


where 


and 


(12) 


(13) 


(14) 
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To  minimize  w  with  respect  to  6,  we  diffe  entiate  w  with  respect  to  y  (since  y  increases  monotonically 
with  <$)  and  set  the  derivative  to  zero.  This  yields  the  following  equation  for  y 


The  solution  to  this  cubic  equation  is 


y3  =  cy  +  1 . 


(15«) 


y  = 


(156) 


which  can  be  evaluated  for  a  given  L  to  yield  the  optimum  6.  Once  y  is  determined  we  can  solve  for  the 
optimum  pixel  spacing,  6a ,  from  the  following  equation: 


K  = 


(16) 


It  can  be  shown  that  f>„  increases  as  L  increases. 

The  above  solution,  however,  may  not  always  be  realizable.  In  any  practical  system,  the  minimum  6 
(denoted  as  6,ntn)  is  limited  by  the  imaging  system  to  a  value  larger  than  the  wavelength  A.  If  we  use  an 
imaging  lens  of  F/number  z,  the  smallest  resolvable  spot  is 


tmin  =  W4-*2  +  1, 

which  corresponds  to  the  highest  spatial  frequency  planewave  traveling  at  an  angle 


(17) 


0,  —  sin-1  (£)• 

Inside  the  crystal,  this  becomes  (from  Snell’s  law) 


(18) 


(19) 


Therefore  the  smallest  resolvable  spot  size  inside  the  crystal  is  also  ^m;n  as  given  by  Eq.  (17). 

f>min  is  the  lower  bound  for  the  size  of  b.  If  L  is  too  small,  S„  (from  Eq.  (15))  becomes  less  then  Amin .  In 
that  case,  we  set  6  =  6m,n.  Since  6„  increases  as  L  increases,  we  can  use  Eqs.  ( 1 3)— ( 16)  to  find  the  smallest 
L  for  which  ba  is  larger  than  6mi„: 


where 


(20) 


(21) 
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\\V  will  refer  to  the  condition  where  the  optimum  :  ss  than  d„„n  as  the  “thin"  disk  regime.  If,  on  the 
other  hand,  the  optimum  pixel  spacing  is  larger  than  the  resolution  limit  of  the  lens,  this  corresponds  to 
the  "thick"  regime.  For  example,  for  n  =  2.2,  A  =  500  inn,  .V,,  =  1000,  and  using  a  F/3  imaging  lens  (i.e., 
;  =  3),  we  get  =  3.04  /mi  and  Ltn,n  =  40. 3G  nun.  Note  that  L„„„  does  not  depend  on  <?•>. 

To  summarize,  if  L  >  Lm,n  (thick  disk),  we  use  <*>  =  />„  (from  Ecj.  (15)),  and  if  L  <  L„lln  (thin  disk),  we 
use  <S  =  (as  given  by  E<i-  (17)). 


2d.  Optimum  Thickness  L 


Our  problem  is  now  reduced  to  maximizing  N  with  respect  to  the  two  remaining  variables  L  and  0j. 
We  first  treat  0 j  as  fixed,  and  find  the  optimum  L  that  maximizes  X/A. 

In  the  range  L  <  Lmtn  we  ttse  A„n„  as  the  optimum  />,  and  write  X/A  from  Eq.  (11)  as 


where 


and 


X/A  = 


1  1  +  itx 

*m.»  (1  +7J-) 


ii  L 

aa;’ 

nX„ 


n  = 


(cos^i  —  cosfr>)> 

1 


/  3 - 

Vmin  V  Vmm 


(22) 


(23) 


(2.) 


(25) 
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tan  02 

i/’ 

«vmm 


(20) 


We  can  solve  for  the  optimum  L  by  differentiating  the  expression  in  Eq.  (22)  with  respect  to  x.  The  maximum 
X/A  turns  out  to  be 


X/A 


o/J-) 


-  (v/W-x/T7!)’ 


(27) 


which  occurs  at 


(28) 


assuming  of  course  that  L„  <  Lmtn.  If  L„  >  Lmi„  this  means  that  the  optimum  thickness  is  out  of  the 
thin  regime  where  the  analysis  used  to  derive  L„  applies.  Within  the  thin  regime  the  maximum  thickness  L 
occurs  at  the  boundary  since  N/A  is  monotonically  increasing  with  L  for  L  <  L,„,n.  To  obtain  the  overall 
optimum  thickness  L  we  must  compare  the  maximum  obtained  from  this  regime  (i.e.,  L  <  L,ntn)  with  the 
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optimum  thickness  obtained  from  the  thick  regime  (£  >  £mi„)  and  finally  select  the  thickness  that  yields 
the  larger  density  N/A. 

As  an  example,  we  continue  with  tne  previous  example  where  (5m;„  =  3.04  //m  (for  a  lens  with  F/number 
of  3).  If  we  take  8 1  =  10°  and  82  =  20°,  we  find  £„  to  be  16.74  mm,  which  is  less  than  £m,„  =  40.36  mm. 
Therefore,  the  solution  obteiu^d  bom  the  thin  regime  is  the  valid  optimum  thickness.  Note  that  as  Np 
increases,  so  does  a,  and  t.Vr-Tore  the  expression  in  Eq.  (27)  increases.  For  large  Np,  the  maximum  N/A  as 
given  by  Eq.  (27)  increase  -  approximately  linearly  with  Np,  or  the  square  root  of  the  total  number  of  pixels 

For  £  >  Lmin ,  we  can  use  Eq.  (16)  for  6  and  using  Eqs.  (1C  1— ( 15),  Eq.  (11)  can  be  written  as 


w-(5)' 


1  +  or3'2 


(29) 


The  above  expression  can  be  evaluated  numerically  to  find  the  value  of  £  which  maximizes  N/A.  We  can 
also  derive  a  relatively  simple  asymptotic  expression  (for  large  £)  by  observing  that  as  £  — >  00,  y  — *  y/c. 
The  asymptotic  expression  for  N/A  is 


...  pi\2  7rA’t>(cos#i  -  cos 82)  I XNP 

1  '  VA/  16 tan #2  V  »£  ' 


(30) 


The  above  expression  predicts  that  the  density  will  decrease  as  the  disk  thickness  becomes  very  large.  This 
is  confirmed  by  the  numerical  results  we  present  in  the  following  section. 


2e.  Optimum  8n  and  the  Maximum  Storage  Density 

The  final  step  in  the  optimization  of  the  storage  density  N/A,  consists  of  optimally  selecting  82 ■  Since 
we  cannot  analytically  derive  the  optimum  angle,  we  resort  to  numerical  methods.  In  Figure  5  we  plot 
Eq.  (22)  in  the  thin  regime  (solid  line)  and  Eq.  (29)  in  the  thick  regime  (dotted  line)  as  a  function  of  £  for 
various  values  of  82  using  the  optimum  value  for  <“>.  The  vertical  line  indicates  the  transition  from  one  regime 
to  the  other.  The  optimum  values  for  £  and  8  are  those  that  yield  the  maximum  density.  The  parameters 
used  in  plotting  Fig.  5  are  A  =  500  nin,  Np  =  1,000,  n  -  2.2  (the  index  of  refraction  for  LiNb03  crystals), 
and  8 1  =  10°.  8 2  =  30°  is  the  maxi.,,  11m  value  for  which  N/A  is  plotted  since  27.04°  is  the  largest  angle 
that  can  be  supported  inside  the  crystal  (due  to  Snell’s  law)  without  resorting  to  the  use  of  index  matching 
fluids. 

From  Figure  5,  we  see  that  the  maximum  N/A  is  obtained  near  £  =  1.5  cm  and  increases  monotonically 
with  82  for  the  parameters  we  selected.  In  this  case,  the  optimum  thickness  is  in  the  thin  regime  (£„  <  £m,„). 
Since  it  is  not  practical  to  use  82  =  30°  inside  the  crystal  (the  critical  angle  is  27.04°),  we  get  a  realistic 
estimate  for  the  achievable  density  by  using  82  —  20°.  The  corresponding  angle  swing  outside  the  crystal 
is  then  22.5°  to  48.8°  (total  angular  swing  of  26.3°)  which  is  practically  achievable.  The  maximum  density 
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N/A  for  various  values  of  ©2 


Figure  5.  Angle  multiplexing:  N/A  vs.  L  for  various  values  of  02.  We  take  Np  =  1000,  n  =  2.2, 

A  =  500  nm,  and  —  10°. 

N/A  is  29.3  bits/pm2,  which  is  obtained  for  a  crystal  thickness  of  L  =  16.74  mm  using  Ng  =  1306  angularly 
multiplexed  holograms.  This  density  can  be  increased  by  a  factor  of  4  (giving  us  N/A  =  117.2  bits/pm2 )  if 
we  simultaneously  record  reflection  and  transmission  holograms  in  the  same  reference  angle  range  from  both 
sides  of  the  signal  beam,  as  shown  in  Figure  3.  The  area  for  each  recording  location  is  w  x  w'  =  4.3  x  10.4  mm2. 
Figure  6  is  a  plot  of  the  optimum  density  and  also  the  number  of  angularly  multiplexed  holograms,  Ng ,  as 
a  function  of  L.  For  the  thickness  that  yields  maximum  density,  Ng  =  1,306  holograms.  Since  more  than 
5,000  holograms  have  been  recorded  and  faithfully  reproduced  in  Lithium  Niobate  [10],  the  geometric  factors 
considered  in  this  report  limit  the  recording  more  severely  than  the  material  dynamic  range.  As  another 
example,  if  we  record  only  100  holograms  at  each  location,  then  the  optimum  thickness  is  a  little  over 
1  mm  and  the  corresponding  storage  density  is  about  8.8  bits/pm2  (compared  to  about  30  bits/pm2  for  the 
optimum  design).  This  density  can  be  increased  by  a  factor  of  4  as  we  already  described  in  Figure  3. 


3.  Wavelength  Multiplexing 

Wavelength  multiplexing  [11-12]  is  an  alternative  method  for  multiplexing  holograms  in  a  single  location 
on  the  HD.  In  this  section  we  calculate  the  capacity  of  a  wavelength  multiplexed  HD  using  a  similar  derivation 
as  for  angular  multiplexing.  The  number  of  bits  that  can  be  stored  is  expressed  as 


N  =  N.NxN},  (31) 

where  N\  is  the  number  of  wavelength  multiplexed  holograms.  We  assume  that  the  wavelength  A  sweeps 
from  A!  to  Aj,  with  Ai  <  A2. 
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iVp=1000,  n=2.2,  A=500  nm 


Figure  8.  Angle  multiplexing:  optimum  N/A  (optimized  with  respect  to  <*>)  and  ,Yp  as  functions 
of  thickness  L.  We  take  Np  =  1000,  n  =  2.2,  A  =  500  nm,  6\  =  10°,  and  02  =  20°. 


Figure  7.  Wavelength  multiplexing:  extra  area  taken  up  because  of  defocusing  and  wavelength 
change. 


For  wavelength  multiplexing,  we  again  assume  that  the  image  is  at  normal  incidence  and  focused  at  the 
middle  of  the  crystal  (Figure  7),  and  that  the  reference  beam  is  counter-propagating  with  the  image  beam 
also  at  normal  incidence.  In  this  case,  the  problem  of  image  defocusing  at  the  crystal  surface  is  the  same, 
and  we  get  Eq.  (8)  for  the  width  w  as  before.  However,  since  the  reference  beam  is  co-linear  with  the  signal 
beam  for  all  wavelengths,  there  is  no  extra  width  taken  up  by  the  L  tan  02  term  in  the  expression  for  w'  in 
Eq.  (9).  On  the  other  hand,  as  A  sweeps  through  Aj  to  A2,  u>  changes.  For  any  choice  of  6,  the  largest  w  is 
for  A  =  A2.  Therefore  we  have 


12 


(32) 


where 

...  c  L 

w  =  A  +  — 7 . . . . 

P  s/{n6/ \2)2-l 

is  a  function  of  Aj. 

To  find  N\,  we  note  that  the  half-width  of  the  (frequency)  selectivity  A u  is  [8] 


(33) 


(34) 

where  t>c  is  the  speed  of  light  in  vacuum.  As  A  sweeps  over  Ai  to  A2,  the  number  of  wavelength  multiplexed 
holograms  that  can  be  stored  is  therefore 


Nx  =  1  +  HjU1  =  1  + 

2Au 


nL  M  _  _1_\ 

2  VAi  >2/’ 


(35) 


where  we  take  the  separation  between  adjacent  holograms  to  be  a  full  width  2Au.  Using  Eqs.  (32),  (33), 
and  (35),  we  then  have 


N  =  AN* 


>  +  t(i-i) 


P  r 


Np6  + 


v/(ni/A,)3-l 


2  • 


(36) 


3a.  Optimum  Np,  Aj,  and  6 

We  now  want  to  maximize  N  with  respect  to  Np,  L,  Aj,  A2,  and  6.  As  before,  N  increases  monotonically 
with  Np,  which  is  limited  by  the  SLM  to  about  1,000.  N  also  increases  as  the  minimum  wavelength  Aj 
decreases.  This  will  be  limited  by  the  shortest  usable  wavelength  we  can  get  out  of  a  tunable  laser  and/or 
the  spectral  sensitivity  of  the  material.  For  the  remainder  of  this  section,  we  will  assume  that  Np  and  Ai  are 
given  and  fixed. 

The  three  remaining  parameters  A2,  L,  and  6  are  more  complicated.  We  first  take  L  and  A2  as  fixed, 
and  find  the  optimum  6.  Considering  N/A  as  a  function  6,  we  find  as  before  that  the  maximum  N/A  is 
obtained  when  w  is  minimized  with  respect  to  6.  We  then  get  the  same  set  of  equations  as  Eqs.  (13)-(16), 
except  with  A  replaced  by  A2.  We  also  have  the  same  6rai„  and  Lmi„  (with  A  replaced  by  A2)  conditions  as 
given  respectively  by  Eqs.  (17)  and  (20).  Note  that  both  6mi„  and  £m,„  scale  linearly  with  wavelength  (since 
!/min  depends  only  on  n  and  z,  the  F/number  of  the  imaging  lens).  It  should  be  emphasized,  that  is 
the  resolution  of  the  imaging  system  using  wavelength  A2.  The  resolution  of  the  system  using  Ai  (which  is 
less  than  A2)  is  of  course  better. 

In  summary,  if  L  >  im;„,  we  use  6  =  6a  (from  Eq.  (15)),  otherwise  we  use  6  =  (as  defined  in 
Eq.  (17));  in  these  equations,  A  is  replaced  by  A2.  As  an  example,  for  A2  =  540  nm  and  a  imaging  lens 


13 


with  F/numoer  of  3,  we  have  6min  =  3.28  pm  and  £m;n  =  43.59  mm.  For  A2  =  750  nm,  these  become 
6min  =  4.56  pm  and  Lmtn  =  60.54  mm. 

3b.  Optimum  L  and  A2 

We  now  find  the  optimum  thickness  L  that  maximizes  N/A.  In  the  thin  regime  ( L  <  Lmtn)  we  take 
6  =  and  write  N/A  as 


where 


N/A  = 


1  1  +  ax 

d(l  +  ^)2’ 


and 


(37) 


(38) 

(39) 


/?  = 


ymin  Vrfn.n-1 

By  differentiating  the  expression  in  Eq.  (37)  with  respect  to  x,  we  find  the  maximum  N/A  to  be 


(40) 


N/A  = 


I 


a 


1 


e,n  4  0(0 -ay 


(41) 


which  occurs  at 

For  example,  for  Aj  =  500  nm,  A2  =  540  nm  (A2/A1  =  1.08),  n  —  2.2,  and  Np  =  1000,  we  get  L0  =  43.82  mm. 
If  A2  increases  to  750  nm,  La  increase  to  60.88  mm.  In  both  cases,  L„  is  larger  than  Lmi„  (43.59  mm  and 
60.54  mm  respectively).  This  means  that  in  there  is  no  maximum  in  the  thin  regime  and  therefore  in  the 
range  L  <  Lmi„,  N/A  is  monotonically  increasing  with  L.  In  this  case,  we  would  select  the  boundary 
value  ( Lmin )  for  the  best  thickness  obtainable  from  the  thin  regime.  Notice  that  for  wavelength  multiplexed 
storage  the  optimum  thickness  of  the  disk  can  become  quite  large.  Even  though  we  are  not  considering 
materials  issues  in  this  report,  we  should  point  out  that  the  useful  thickness  of  the  material  in  practice  can 
be  limited  by  absorption.  In  some  materials  (e.g.,  Lithium  Niobate)  it  is  possible  to  reduce  the  absorption  by 
properly  preparing  the  material  (e.g.,  by  adjusting  the  dopant  and  reduction/oxidation  level).  The  reduced 
absorption  will  typically  reduce  the  recording  speed  of  the  material  for  a  given  light  intensity.  Therefore, 
when  materials  considerations  are  included  in  the  design  process,  this  trade-off  between  speed  and  density 
will  emerge. 
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In  the  thick  regime,  L  >  Lmi„,  we  use  Eqs.  (15)  and  (16)  (with  A  replaced  by  A2)  to  obtain  6,  and  write 
N/A  as 


N/A  = 


1  -f  ac3! 2 


(43) 


where  a  is  given  by  Eq.  (39).  As  before,  as  L  — *  00,  y  — >  y/c.  For  wavelength  multiplexing,  however,  the 
asymptotic  behavior  of  N/A  is  different.  As  L  — *  00,  N/A  saturates  and  approaches 


N/A 


n2N„ 


8  A 


W  \Ai  A2/ 


(44) 


Thus  for  the  range  L  >  Lmi„ ,  N/A  also  increases  monotonically  with  L. 

Note  that  the  saturation  value  of  N/A  increases  as  Np  increases  and  Ai  decreases.  Also,  for  any  choice 
of  A2/At,  we  have 


±(±.±)<±. 
a2  U  a2;-4a r 


with  equality  at 


(45) 


(46) 


Thus,  even  if  it  is  practical  to  have  a  light  source  with  such  a  large  range  of  wavelength  tunability,  the 
optimum  setting  for  A2/Ai  in  order  to  obtain  maximum  saturation  density  is  2  (provided  we  use  the  same 
Ai).  For  practical  systems,  A2/At  is  smaller  than  2,  and  in  this  range  the  saturation  value  of  N/A  increases 
as  A2/Ai  increases.  In  the  case  where  AA  =  A2  -  Ai  <  Ai  (A2  =s  Ai),  the  saturation  value  given  in  Eq.  (44) 
is  approximately 


N/A  x 


n2Np  AA 
8  A  ?’ 


(47) 


which  is  proportional  to  AA. 

In  practice,  the  range  of  usable  wavelengths  is  determined  by  the  laser  system.  For  instance,  dye  lasers 
can  be  tuned  in  the  range  from  370  nm  to  890  nm,  which  gives  us  a  A2/Ai  of  2.40,  in  excess  of  the  optimum 
A2/A1  =  2  requirement.  It  should  be  noted,  however,  that  it  is  necessary  to  use  several  different  dyes  in 
order  to  achieve  this  range  of  wavelengths.  For  a  typical  broadband  laser  dye  such  as  Coumarin  6,  the  range 
is  from  510  nm  to  550  nm,  which  only  gives  us  a  A2/Aj  of  1.08.  For  TirSapphire  lasers,  the  range  is  690  nm 
to  1025  nm,  which  gives  us  a  A2/Ai  of  1.48. 

As  a  specific  example,  consider  the  case  where  Np  =  1000,  »  =  2.2,  and  At  =  500  nm.  We  plot  N/A  as 
a  function  of  L  (where  N/A  has  been  optimized  with  respect  to  <“>)  for  various  values  of  A2/Aj  .  The  result  is 
shown  in  Figure  8.  We  see  that  N/A  saturates  for  large  L  (around  5  cm)  as  expected,  and  the  saturation 
value  is  largest  for  A2/Ai  =  2.  In  Figure  9,  we  plot  N/A  and  N\  as  functions  of  L  for  A2/Ai  =  1.08  and 
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N/A  for  various  values  of  X2/Xj 


L  (mm) 

Figure  8.  Wavelength  multiplexing:  optimum  N/A  (optimized  with  respect  to  A)  as  a  function  of 
L  for  various  values  of  A2/Ai.  We  take  At  =  500  nm,  Np  =  1000,  and  n  =  2.2. 


.  A  =1000,  n= 2.2.  X,=500  nm 
10s  p 

104 

10s 

102 

10' 

10° 

10'2  10'1  10°  10'  10*  103 
L  (mm) 

Figure  0.  Wavelength  multiplexing:  optimum  N/A  and  N\  as  functions  of  thickness  L  for  A2 /X\  — 

1.08  and  A2/Aj  =  1.5.  We  take  Ai  =  500  nm,  Np  =  1000,  n  =  2.2. 

A2 /Aj  =  1.5  using  the  same  Np,  n,  and  Aj.  For  A2/Ai  =  1.08,  N/A  approaches  166.0  bits/pm2,  while  for 
A2 /A j  =  1.5,  N/A  approaches  537.8  bits//rm2. 

3c.  Storage  Density  and  Optimum  A2  for  “Thin”  Disks 

The  point  where  L  causes  N/A  to  reach  saturation  is  of  the  order  of  5  cm.  At  this  thickness,  it  becomes 
questionable  about  what  we  mean  by  a  “disk”.  In  practice  it  may  be  desirable  or  necessary  (e.g.,  because  of 


16 


absorption)  to  keep  the  thickness  small.  In  this  case  we  are  in  the  L  <  Lmi„  range  (even  though  L0  may  be 
larger  than  Lm,„),  and  N/A  is  given  by  Eq.  (37).  We  can  approximate  Eq.  (37)  by 


if  we  assume  that 


(48) 


ax  >  1  >  (3x .  (49) 

In  the  previous  example,  a  =  40  for  A2/A1  =  1.08,  and  a  =  250  A2/A1  =  1.5,  while  (3  =  5.60  x  10~3  in  both 
cases,  so  the  condition  is  satisfied. 

If  we  limit  x  and  hence  the  disk  thickness  L  to  the  range  required  by  Eq.  (49),  the  optimum  A2  can  be 
found  by  taking  the  derivative  of  the  expression  in  Eq.  (48)  with  respect  to  Ao  and  setting  it  to  zero.  In  this 
case,  it  is  easy  to  show  that  the  maximum  N/A  occurs  for  A2/A1  =  1.5  (again  assuming  that  we  are  using 
the  same  At),  which  is  very  close  to  the  value  provided  by  Ti:Sapphire  lasers.  Therefore,  in  this  case  the 
density  does  not  increase  indefinitely  with  A  A. 

Finally  we  can  also  calculate  the  “knee”  of  the  N/A  curve,  which  we  define  as  the  point  were  the 
expression  given  by  Eq.  (47)  reaches  the  saturation  value.  This  is  given  by 


L  =  Lk 


n  NPSjtn 
4A2 


4r2  -f  1 

4 


llX^Ny, 


(50) 


which  is  proportional  to  A2.  For  A2 /A  1  =  1.08,  Lk  =  11.0  mm,  which  gives  us  N/A  =  106.5  bits//im2  and 
N\  =  1,794.  For  A2/A1  =  1.5,  Lk  —  15.3  mm,  which  gives  us  N/A  =  34.46  bits//rm2  and  N\  =  11,221.  In 
both  cases,  Lk  is  less  than  iml„,  and  the  corresponding  values  of  N/A  is  slightly  over  half  of  the  saturation 
values  for  N/A  (i.e.,  approximately  a  3  dB  drop). 


4.  Discussion  And  Conclusions 

The  values  for  the  various  parameters  discussed  in  this  report  are  summarized  in  Table  I,  and  the  storage 
densities  N/A  are  plotted  in  Figure  10  where  we  denote  the  densities  of  angle  multiplexing  and  wavelength 
multiplexing  by  ( N/A)g  and  (N/A)\,  respectively. 

In  Figure  10,  the  curves  for  (N/A)g  using  just  the  angle  range  Q\  to  O2  (either  as  transmission  or 
reflection  holograms)  is  marked  as  ( x  1).  We  see  that  it  is  about  a  factor  of  2  smaller  than  (N/A) However, 
if  we  angle-multiplex  from  both  sides  of  the  signal  beam,  ( N/A)»  increases  by  a  factor  of  2  (denoted  by  the 
x2  curve  in  Figure  10).  If  we  further  record  both  reflection  and  transmission  holograms  (as  in  Figure  3), 
this  increases  by  a  factor  of  4  (denoted  by  the  x4  curve  in  Figure  10).  In  this  case,  (N/A)g  becomes  larger 
than  (N/A)x  until  L  reaches  about  12.5  mm,  where  both  (N/A)\  and  (N/A)g  are  about  115  bits//im2. 
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TABLE  I.  Value  of  Parameters  used  in  Figure  10 


Parameters 

Index  of  Refraction 
Number  of  Pixels 
Wavelength 
Angles 
Pixel  Size 
Critical  Thickness 
Optimum  Thickness 
Maximum  Density 
Number  of  Holograms 


(F /number  of  imaging  lens  is  3) 

Angle  multiplexing  Wavelength  multiplexing 


n  =  2.2 
JV2  =  106 
A  =  500  nm 

0i  =  10°,  0!  =  20°, 

Smin  =  3.04  pm 

Lmin  =  40.36  mm 

L„  =  16.74  mm 

N/A  =  4x  29.3  bits/pm2 

Ne  =  1306 


n  =  2.2 

V2  =  106 

A]  =  500  nm,  A2  =  540  nm 

Smi„  =  3.28  pm 
Lmin  =  43.59  mm 
La  ss  30  mm 
N/A  =  166.0  bits/pm2 
N\  ~  5000 


L  (mm) 

Figure  10.  Comparison  of  angle  multiplexing  and  wavelength  multiplexing.  We  take  A  =  Aj  = 
500  nm,  A2/A1  =  1.08,  n  =  2.2,  Np  =  1000,  0i  =  10®,  and  02  =  20°.  The  density  using  angle 
multiplexing  is  denoted  by  ( N/A)g ,  and  the  density  using  wavelength  multiplexing  is  denoted  by 
(N/A)  A. 


One  might  ask  whether  it  is  possible  to  achieve  higher  density  by  recording  in  the  Fourier  plane  instead 
of  the  image  plane.  It  turns  out  that  the  storage  density  is  the  same.  This  is  because  the  space- bandwidth- 
product  is  a  constant.  Specifically,  consider  an  image  of  extent  a  =  Np6,  where  6  is  the  pixel  spacing  and 
Np  is  the  number  of  pixels  along  one  dimension.  Let  b  be  the  extent  of  the  Fourier  transform  of  this  image 
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(by  a  lens  of  focal  length  F),  and  let  1/S'  be  the  highest  spatial  frequency  of  the  Fourier  transform.  Then 
within  the  paraxial  approximations 


Np  = 


(51) 


This  shows  that  recording  in  either  the  image  plane  or  the  Fourier  plane  will  give  the  same  minimum  width. 

If  we  record  holograms  at  off-image  or  off-Fourier  planes,  the  required  width  w  increases.  However,  it 
is  sometimes  desirable  to  do  this  for  purpose  of  noise,  image  quality,  and  alignment  sensitivity.  The  tradeoff 
between  these  requirements  and  storage  density  will  have  to  be  considered  in  the  design  of  a  practical  system. 

We  have  derived  the  optimum  conditions  for  achieving  the  maximum  storage  density  of  a  3-D  HD 
disk  using  either  angle  multiplexing  or  wavelength  multiplexing.  Such  optimally  designed  disks  can  store 
information  with  area  densities  more  than  100  bits/pm2  with  disk  thickness  approximately  1  mm.  However, 
the  limits  to  storage  density  derived  in  this  report  are  only  due  to  the  geometry  of  the  system.  The  storage 
density  can  also  be  limited  by  noise  (cross-talk,  detector  noise,  media  defects,  etc.)  and  the  limited  dynamic 
range  of  the  recording  medium.  These  limits  to  N/A  (addressed  in  a  separate  report  [13])  prove  less  restrictive 
than  the  geometric  limits  derived  here.  This  is  supported  by  recent  experiments  [14]  where  1000  holograms 
were  superimposed  and  reconstructed  with  extremely  low  probability  of  error  in  a  lithium  niobate  crystal 
with  1  cm  thickness.  The  parameters  of  this  experiment  were  reasonably  close  to  the  optimum  parameters 
we  derived. 


19 


REFERENCES 


[1]  P.J.  Van  Heerden,  “Theory  of  Optical  Information  Storage  in  Solids,”  Appl.  Optics ,  2,  393-400  (1963). 

[2]  K.  Blptekjaer,  “Limitations  on  Holographic  Storage  Capacity  of  Photochromic  and  Photorefractive 
Media,”  Appl.  Opt.  18,  57-67  (1979). 

[3]  D.  von  der  Linde  and  A.M.  Glass,  “Photorefractive  Effects  for  Reversible  Holographic  Storage  of  Infor¬ 
mation,”  Appl.  Phys.,  8,  85-100  (1975). 

[4]  J.-P.  Huignard,  J.-P.  Herriau,  and  F..  Micheron,  “Coherent  Selective  Erasure  of  Superimposed  Volume 
Holograms  in  LiNbOa,”  Appl.  Phys.  Lett.,  26,  256-258  (1975). 

[5]  L.  d’Auria,  J.-P.  Huignard,  C.  Slezak,  and  E.  Spitz,  “Experimental  Holographic  Read-Write  Memory 
using  3-D  Storage,”  Appl.  Opt.,  13,  808-818  (1974). 

[6]  Y.  Nakagome,  H.  Tanaka,  K.  Takeuchi,  E.  Kume,  Y.  Watanabe,  T.  Kaga,  Y.  Kawamoto,  F.  Murai, 
R.  Izawa,  D.  Hisamoto,  T.  Kisu,  T.  Nishida,  E.  Takeda,  and  K.  Itoh,  “An  Experimental  1.5-V  64-Mb 
DRAM,”  IEEE  J.  Solid  State  Circuits,  26(4),  465-472  (1991). 

[7]  D.  Psaltis, “Parallel  optical  memories”,  Byte,  17(9),  179-182  (1992). 

[8]  H.  Kogelnik,  “Coupled  Wave  Theory  for  Thick  Hologram  Gratings,”  Bell  Syst.  Tech.  J.,  48(9),  2909- 
2947  (1969). 

[9]  C.  Gu,  J.  Hong,  I.  McMichael,  Ragini  Saxena,  and  F.H.  Mok,  “Cross  Talk  Limited  Storage  Capacity  of 
Volume  Holographic  Memory,”  J.  OS  A.  A,  9,  1-6  (1992). 

[10]  F.H.  Mok,  “Applications  of  Holographic  Storage  in  Lithium  Niobate,”  (presented  at  OSA  1992  Annual 
Meeting)  OSA  1992  Annual  Meeting  Technical  Digest,  Vol.  23,  WEI,  102,  Sept  1992. 

[11]  F.T.S.  Yu,  S.D.  Wu,  A.W.  Mayers,  and  S.M.  Rajan,  “Wavelength  Multiplexed  Reflection  Matched 
Spatial  Filters  Using  LiNbOa,”  Opt.  Comm.,  81(6),  343-347  (1991). 

[12]  G.A.  Rakuljic,  V.  Leyva,  and  A.  Yariv,  “Optical-data  Storage  by  using  Orthogonal  Wavelength- 
Multiplexed  Volume  Holograms,”  Opt.  Lett.,  17(20),  1471-1473  (1992). 

[13]  H.-Y.  Li  and  D.  Psaltis,  “3-D  Holographic  Disks:  Material  Considerations,”  to  be  published. 

[14]  F.  Mok,  D.  Psaltis,  and  G.  Burr,  “Spatial  and  Angle  Multiplexed  Holographic  Random  Access  Memory,” 
SPIE  vol.  1773c,  San  Diego,  CA,  July  1992. 


20 


FIGURE  CAPTIONS 


1.  The  3-D  holographic  disk  (HD). 

2.  Recording  geometry. 

3.  Angular  multiplexing  by  reflection  and  transmission  holograms  from  both  sides  of  the  signal  beam. 

4.  Angle  multiplexing:  extra  area  taken  up  by  defocusing  and  reference  beam  angle  change. 

5.  Angle  multiplexing:  N/A  vs.  L  for  various  values  of  We  take  Np  =  1000,  n  =  2.2,  A  =  500  nm,  and 
9X  =  10°. 

6.  Angle  multiplexing:  optimum  N/A  (optimized  with  respect  to  6)  and  Np  as  functions  of  thickness  L. 
We  take  Np  =  1000,  n  =  2.2,  A  =  500  nm,  9X  =  10°,  and  9 2  =  20°. 

7.  Wavelength  multiplexing:  extra  area  taken  up  because  of  defocusing  and  wavelength  change. 

8.  Wavelength  multiplexing:  optimum  N/A  (optimized  with  respect  to  6)  as  a  function  of  L  for  various 
values  of  A2/Ai.  We  take  At  =  500  nm,  Np  —  1000,  and  n  =  2.2. 

9.  Wavelength  multiplexing:  optimum  N/A  and  N\  as  functions  of  thickness  L  for  \-i/\x  —  1.08  and 
A2/Ai  =  1.5.  We  take  Ai  =  500  nm,  Np  =  1000,  n  =  2.2. 

10.  Comparison  of  angle  multiplexing  and  wavelength  multiplexing.  We  take  A  =  Aj  =  500  nm,  A2/Aj  = 
1.08,  n  =  2.2,  Np  —  1000,  9X  =  10°,  and  92  —  20°.  The  density  using  angle  multiplexing  is  denoted  by 
(N/A)g,  and  the  density  using  wavelength  multiplexing  is  denoted  by  {N / A) x- 
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